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Abstract. In this review, we highlight the physics of synchronization 
in collections of beating cilia and flagella. We survey the theory syn¬ 
chronization in collections of noisy oscillators. This framework is ap¬ 
plied to flagellar synchronization by hydrodynamic interactions. The 
time-reversibility of hydrodynamics at low Reynolds numbers prompts 
swimming strokes that break symmetry to facilitate hydrodynamic syn¬ 
chronization. We discuss different physical mechanisms for flagellar syn¬ 
chronization, which break this symmetry in different ways. 


1 Introduction 


What we call today synchronization of oscillators was first described in 1665 by van 
Huygens, when he observed two pendulum clocks adopting a common rhythm that 
lasted for many hours . This was surprising, given the notorious inaccuracy of clocks 
at the time. In fact, a weak mechanical coupling had entrained the two pendulum 
clocks. A simple class-room demonstration of this experiment can be realized using 
two metronomes on a swinging tray [^. In the middle of the 20*^ century, work on 
phase-locked electronic oscillators motivated the development of a general theory of 
synchronization |^|^, which we briefly review in section 2. 

Synchronization of active oscillators has since then been observed in a number of 
examples, including synchronization of biological oscillators. As a popular example, 
synchronization of walking gaits among pedestrians was observed on the Millennium 
bridge in central London, which caused large-amplitude vibrations of the bridge mUm . 
This synchronization was the result of a mechanical coupling. Walking of individual 
pedestrians excited small-amplitude vibrations; in turn, each pedestrian adapted its 
frequency and phase of walking to the vibrations of the bridge. 

Synchronization of biological oscillators is also observed at the cellular scale. Many 
cells are propelled in a liquid by regular bending waves of slender cell appendages 
termed cilia and flagella [7}|^. The rhythmic beat of a cilium or flagellum represents 
a prime example of a biological oscillator. Cilia and flagella share a highly conserved 
core, the axoneme, which comprises a regular cylindrical arrangement of 9 double t 
microtubules, which are connected by 10^ — 10® dynein molecular motors 
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These molecular motors convert chemical energy in the form of ATP into mechanical 
work to bend the axoneme 
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A dynamic instability of collective motor dynamics 
gives rise to self-organized beiiSing waves of cilia and flagella [T3}|T5) . The distinction 
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between eukaryotic cilia and eukaryotic flagella is mainly made with respect to their 
length and function, and we will use the same term flagellum for both. Note, however, 
that the eukaryotic flagellum is very different from the bacterial flagellum, which is a 
passive protein filament. 

Collections of flagella can adopt a common rhythm and synchronize their beat. 
It had been observed already by Gray almost 100 years ago that pairs of sperm cells 
swimming in close proximity can synchronize their flagellar beat 16,17. At high 


sperm densities, swimming sperm can self-organize into vortex patterns with vortices 
of about 5 cells with mutually phase-locked flagellar beat Flagellar synchroniza¬ 
tion plays an important role for the collective dynamics in carpets of short flagella 
termed cilia, which can phase-lock their beat to give rise to propagating metachronal 
waves. This coordinated flagellar beating results in effective surface slip and thereby 
facilitates self-propulsion of multi-ciliated organisms such as the unicellular eukaryote 
Paramecium 19 or the green alga colony Volvox 20 . Inside the mammalian body. 


dense carpets of cilia on epithelial surfaces p um p fluids, such as mucus in mammalian 
airways or cerebrospinal fluid in the brain 21 . Again, these cilia synchronize their 


beat to exhibit metachronal waves. This self-organized ciliar dynamics is important 
for efficient fluid pumping 22 ■ 24 . 


The bi-flagellated, unicellular green alga Chlamydomonas is serving as a model 
organism to study physical mechanisms of flagellar synchronization 25 , which is 


singled out by the ease of its experimental handling and accessibility for quantitative 
study 26 -30 . Chlamydomonas swims like a breast-stroke swimmer with two flagella, 


which can synchronize their beat, both in free-swimming cells and cells held in a micro¬ 
pipette. On physical grounds, this flagellar coordination is surprising, since there is no 
direct coupling between the respective molecular motors that drive the beat of each 
flagellum. There is also no evidence for a chemical master oscillator that could set a 
common rhythm for both flagella. Below, we discuss different physical mechanisms 
that can account for the stabilization of synchronous flagellar beating by a mechanical 
coupling. In a generic description of synchronization in pairs of oscillators known as 
the Adler equation discussed below, there are two possible synchronization states: in- 
phase and anti-phase synchronization. In Chlamydomonas^ in-phase synchronization 
corresponds to a mirror-symmetric, breaststroke-like dynamics of the two flagella and 
anti-phase synchronization to a phase-shift of approximately tt between both flagella. 
While wild-type Chlamydomonas cells usually display in-phase flagellar synchrony, 
stable anti-phase synchronization has been observed in a flagellar mutant 31 . For 
Chlamydomonas^ synchronized beating is required to swim both fast and straight 


The flagellar beat responds to external mechanical forces, such as hydrodynamic 
friction forces. It was observed that the beat of sperm cells slowed down upon an 
increase in the viscosity of the swimming medium 32 33 . Friedrich et al. measured a 
dynamic force-velocity relationship of the flagellar beat for Chlamydomonas cells by 
correlating changes of the flagellar phase speed during a beat cycle with a rotational 
motion of these cells, which imparted known hydrodynamic forces on each flagellum 
. Transient external flows were shown to reversibly perturb a limit cycle of flagellar 
bending waves 34 . Using strong external flows it was even possible to stop and 


re-start flagellar oscillations in a reversible manner 35 . This susceptibility of the 


flagellar beat to external forces allows for the synchronization of flagellar beating 
to an oscillatory driving. Flagellar beating could be entrained to periodic mechanical 
forcing using vibrating micro-needles [36| or oscillatory external flows [37] . It has been 
proposed already by Taylor that a weak mechanical coupling betweenflagella could 
also underly the striking phenomenon of spontaneous synchronization in collections 
of cilia and flagella. 
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Taylor considered direct hydrodynamic interactions be twe en flagella as one pos¬ 
sibility for a mechanical coupling between beating flagella 38 . Direct hydrodynamic 
interactions between flagella refer to hydrodynamic flows generated by one flagellum 
that propagate to a nearby flagellum and exert hydrodynamic friction forces on it. 
This concept of hydrodynamic synchronization sparked a continuous research effort, 
both experimental and theoretical. A clear demonstration of synchronization by di¬ 
rect hydrodynamic interactions was achieved only recently in an experimental setup of 
two flagellated cells held by separate micro-pipettes at a distance 39 . In addition to 


synchronization by direction hydrodynamic interactions, a mechanism of mechanical 


self-stabilization was proposed for free-swimming cells 30 40,41 


Important insights into physical mechanisms of hydrodynamic synchronization 
have been gained by the study of artificial micro-swimmers and micro-actuators, such 
as colloids driven by oscillating magnetic fields or ‘light-mill’ micro-rotors driven by 


laser light, which likewise perform cyclic motions 42 


46 . Several of these oscillators 


ing mediated by hydrodynamic 


can synchronize their oscillations through a weak coup! 
interactions. 

Synchronization by hydrodynamic interactions is a subtle phenomenon that re¬ 
quires broken symmetries [^ [47}|^ . The physical reason for this is the time-reversal 
symmetry of the Stokes equation, which governs hydrodynamics in the low Reynolds 
number regime of flagellar flows [^ . As a consequence, active oscillators such as beat¬ 
ing flagella must break either spatial or temporal symmetries for synchronization 
to occur. We will discuss different physical mechanisms for flagellar synchroniza- 
which break this symmetry in different ways. 


47 48 51,52 


tion 

In section 2 of this review, we will discuss the nonlinear dynamics of synchro¬ 
nization of noisy oscillators. In section 3, we will apply this theory in the context of 
flagellar synchronization, and discuss the role of broken symmetries for synchroniza¬ 
tion. 


2 Nonlinear dynamics of synchronization 

2.1 Active oscillators 

We consider a prototype of an active oscillator, the normal form of a Hopf oscillator 
with complex oscillator variable Z = Aexp(i(/?) [53| 

Z = iujZ + ^i{Al-\Z\'^)Z. (1) 

The onset of flagellar oscillations had been previously described as a supercritical Hopf 
bifurcation [54| with normal form given by Equation 0. Note that the nonlinear 
friction term /^(Aq — \Z\'^) becomes negative for jZj < Aq. Negative friction is a 
common feature of active systems characterized by positive feedback loops. In these 
systems, a reservoir of energy becomes depleted to drive the oscillations, with energy 
eventually dissipated as heat. 

At steady-state. Equation Q implies sustained oscillations Z = Aq exp{iip) with 
amplitude Aq. The dynamics of the phase ip is given by 

If = LU. ( 2 ) 

This simple description of a phase oscillator generalizes to a large class of nonlinear 
oscillators, and proves to be particularly useful if the amplitude of oscillations stays 
approximately constant. It should be emphasized that ip is only defined up to multiples 
of 27r and that any physical quantity that depends on ip must in fact be a 27r-periodic 
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Fig. 1. Phase-space of a limit cycle oscillator. A. The red trajectory represents a limit 
cycle, i.e. a stable, periodic orbit of the dynamics. Example trajectories (green) that start in 
the vicinity of this limit cycle will eventually converge onto this limit cycle. The black line 
indicates a so-called isochrone: trajectories that start on the same isochrone will converge 
on the same point on the limit cycle. The flow-held is indicated by gray arrows. B. Same 
as panel A, but for a non-isochronous limit cycle oscillator with non-radial isochrones. The 
simulated dynamics corresponds to a non-isochronous Hopf oscillator with Z = iu)Z -|- (/r -f 
itJi)(Ao — |^p)Z, where the isynchrony parameter cui = 0 in panel A and wi > 0 in panel 
B. C. Noisy limit cycle oscillations of a complex oscillator variable Z, see also Equation (|^. 


function. In real applications, one would compute Z and its phase using either the 
Hilbert transform of x (analytical signal), or by limit-cycle reconstruction as discussed 
now. 

We define an oscillator as a stable limit cycle, i.e. a periodic orbit such that 
x(t -\-T) = x(t) for some period T, for which all trajectories starting in a sufficiently 
thin tubular neighborhood of the orbit will eventually converge to this orbit. This 
is illustrated in Figure Note that every limit cycle can be parameterized by a 
phase variable (/?, such that Lp advances by 2t: during one cycle. We require that the 
phase variable advances with uniform speed along the limit cycle. Given an arbitrary 
parametrization of a limit cycle by a cyclic variable 9, this can always be achieved by 
an appropriate reparametrization ip = p{0) [55[|56| . The phase parameterization is 
unique up to the choice of start point. The phase parameterization can be extended 
to a tubular neighborhood of the limit cycle (by defining so-called iso-chrones) [^. 
Thereby, the nonlinear dynamics in the vicinity of a limit cycle can be mapped on 
the simple phase oscillator description given in Equation ([^. 

This definition generalizes in a straight-forward manner to stochastic dynamical 
systems, where the presence of noise implies that trajectories fluctuate around the 
limit cycle. The variance of amplitude fluctuations is set by a competition of the 
‘attractor strength’ of the limit cycle (Lyapunov exponent) and the noise strength. 
An example of such noisy oscillations is given by a (special case) of a noisy Hopf 
oscillator with complex oscillator variable Z [4|[^ 

Z = iujZ + p{Al-\Z\^)Z + Z{U+i^^). (3) 


where and are uncorrelated Gaussian white noise variables satisfying 
2DjSjkS(t — t'), i,j € {A, (/?}, with phase and amplitude noise strength and Da. 
Here, Stratonovich calculus is used. In the limit of weak noise, we End effective Equa¬ 
tions for amplitude and phase. Amplitude fluctuations are governed by a noisy relax¬ 
ation process that always reverts to the mean amplitude Aq |58| . In the limit of weak 


noise, we obtain an Ornstein-Uhlenbeck process reverting to the mean amplitude Aq, 

A = 2/iAQ(Ao — A) -k Aq^^, (4) 


while the noisy phase dynamics is characterized by phase diffusion with phase diffusion 
coefficient D = 


(fi = oj 


( 5 ) 
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For simplicity, Equation ^ refers to an isochronous oscillator with isochrony pa¬ 
rameter wi = 0. In the general case of a non-isochronous oscillator, the coefficient 
p, has to be replaced by ^ -I- iuji. In this case, an effective phase diffusion coeffi¬ 
cient D = Dip + {ujiAq /is found which combines the contribution from phase 
fluctuations and amplitude fluctuations. 

In applications, the phase diffusion coefficient D can be inferred from the width 
Auj = 2D at half-maximum of the Fourier peak of the power spectral density of Z, or 
from the inverse exponential decay time of the phase correlation function |(exp[f(/j(t-|- 
At) — ii^{t)\)\ = exp{—D At). 


2.2 The Adler Equation for a pair of coupled oscillators 


We consider two phase oscillators, which are coupled by a generic coupling term C{6) 

(pL=i^L + C{(pL-‘PR), ( 6 ) 

‘pR = — (Pl)- (7) 


For mechanical oscillators, such a coupling could arise for example by elastic or hy¬ 
drodynamic interactions. We are interested in the dynamics of the phase-difference 
5 = ipL — Pr- To be physically sound, the coupling term C{5) must be a 27r-periodic 
function of 5. Thus, we can expand C(5) as a Fourier series, C(5) = C'^ cos{n6) + 

C" sin(n(5). We find that the dynamics of 6 depends only on the odd coupling terms, 
S = Auj + 2 C" sin(n(5), where Auj = ujr — ujr. In many applications, C{S) is 
dominated by its principal Fourier mode, allowing us to write 


(5 = Z\a; — Asin(i5). (8) 

Equation is the well-known Adler equation I^, where we have set A = —2(7". In 
the special case of zero frequency mismatch, Aui = 0, we can readily read off the 
steady states of the Adler Equation (|^ as <5 = 0 and S = n. In-phase synchronization 
with (5 = 0 is stable for A > 0 and unstable for A < 0, see also Figure A small 
frequency mismatch |Z\a;| > 0 between both oscillators will detune the steady state 
phase difference to a value S = sin“^(Z\a;/A), provided |Z\cij| < A. If the phase- 
difference becomes too large, \Auj\ > A, Equation Q does not possess a steady state 
anymore. This implies that synchronization cannot occur, corresponding to a case of 
phase drift, where 5 increases without bound. 

We can account for synchronization in the presence of noise by considering the 
stochastic Adler equation 

^ — A sin((5)-I-(9) 


Here, ^(t) is a Gaussian white noise variable satisfying = 2D — where 

D = Di + D 2 denotes the sum of the phase-diffusion coefficients of the individual 
phase oscillators pL and pR. In the presence of noise, the phase-difference S has 
stationary probability distribution p(S) of finite width; for zero frequency mismatch, 
Auj = 0, this distribution reads 

where Iq in the normalization factor denotes a modified Bessel function of the first 
kind. Noise induces occasional phase-slips, during which one of the two oscillators 
will perform an extra cycle. The presence of a frequency mismatch introduces a bias 
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Fig. 2. Dynamic regimes of the Adler equation. The Adler equation for the phase difference 
5 between two coupled oscillators can be interpreted as the dynamics of an over-damped 
particle moving in a tilted periodic potential U{5) = —AujS — A cos 5 [^. A. For positive 
synchronization strength A > 0 and zero frequency mismatch Alo — 0, in-phase synchro¬ 
nization with 5 = 0 is stable, while anti-phase synchronization with 5 = tt is unstable. B. 
For negative synchronization strength A < 0, the stability behavior is reversed: in-phase 
synchronization is unstable, while anti-phase synchronization is stable. C. In the presence of 
noise, the phase difference <5 will fluctuate around its mean position. Occasionally, phase slips 
occur during which the phase difference changes by ±2?!. D. In case the frequency mismatch 
exceeds the synchronization strength, Auj > |A|, no synchronized states can exist, and the 
system displays phase drift with a phase difference that increases (decreases) monotonically. 


in the rates G+ and G_ of phase-slips in the directions S —>■ S 2Tr and 6 —>■ S — 27r, 
respectively 


G± = ^ \hA../D{\lD)\-^e^V (±^) . 


( 11 ) 


Here, lix denotes the modified Bessel function of the first kind with imaginary index. 
Equation (11) implies that the rate of phase-slips is exponentially suppressed for 
G ^ |A| while for G ^ |A| the dynamics of 5 is virtually indistinguishable from phase 
diffusion. For Auj — 0, we note the asymptotic formulas 


G+ = G_ = 


exp(— 2 A/G)/( 27 r) 

DH2^f 


G< |A| 

G> iAr 


( 12 ) 


Many additional analytical results are known, see e.g. the book by Stratonovich [^. 

In conclusion, a generic anti-symmetric coupling between two oscillators implies 
their mutual synchronization, provided this coupling is sufficiently strong to overcome 
both a possible frequency mismatch and the effect of noise. In the minimal model 
given by Equation (|^, there exists an in-phase and an anti-phase synchronized state. 
Which of the two is stable (and thus will be observed in experiments) depends on the 
sign of the synchronization parameter A and is thus a non-generic system’s property 
that depends on the particularities of the given system. 
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2.3 The Kuramoto model of N coupled oscillators 


We now consider the general case of N oscillators with mutual coupling. In this case, 
the Adler Equation ® generalizes to the famous Kuramoto model 


A ^ 

= w* - — ^ sin((^i - ifj), 
^ k^l 


(13) 


with coupling parameter A > 0. Here, we assumed an all-to-all coupling of the oscil¬ 
lators. The Kuramoto model provides a minimal description of oscillator arrays with 
long-range coupling that is analytically tractable. 

We assign a complex oscillator variable Zk = exp(z(^fe) to each phase oscillator 
ipk- The mean of all specihes an order parameter r = \Z\ of phase coherence 


N 


Z = = r exp{iip). 




With this definition, Equation (13) becomes 


(Pi=uji- Arsin(V^ - ipi), 


(14) 


(15) 


i.e. we can replace the coupling between the individual oscillators by a coupling to a 
mean-held. Many analytical results are known, including a complete set of integrals 
of motion for the special case of identical oscillator frequencies 60 . 

We consider the thermodynamic limit oi N —> oo, with a given distribution of 
oscillator frequencies p{uj). Any oscillator Z for which \Auj\ < Xr will phase-lock to 
Z with a constant phase-lag Aip = sin“^[Z\a;/(Ar)]. Here, Auj = w — wq denotes the 
mismatch between the intrinsic oscillator frequency uj and the global synchronization 
frequency ujq = tp. Thus, the stationary probability density of the oscillator phase (p, 
conditioned by the global phase ip, simply reads ps^plip) = S{(p — ip — Aip) in this 
case. Those oscillators witl^Z\a;| > Xr will not phase-lock to Z and display phase drift 
instead. From Equation ([^), one hnds for the stationary probability density of their 
phase Puipl'ip) ^ — ujq = Auj + Xr sin(Z\(p). By inserting the combined stationary 

probability density p{(p\'ip) into Equation ( |14[ ), one obtains an implicit equation for 
the order parameter r that depends on p{uj). For the special case of a Lorentzian 
distribution of oscillator frequencies with location parameter ujq and half-width at 
maximum 7 given by 


p{uj) = - 


1 


TT 72 -I- (w — ujoy 
one can compute the order parameter r explicitly 


(16) 


A < Ac 

^1- A>Ac. 


(17) 


Here, Ac denotes a critical coupling strength given by Ac = 27 . We conclude that the 
system of coupled oscillations exhibits a phase-transition as a function of synchro¬ 
nization strength A, characterized by a complete lack of synchronization below the 
critical coupling strength, and partial synchronization above. 
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3 Physical principles of flagellar synchronization 


We now apply the generic concepts introduced above to the specific case of beating 
cilia and flagella, which can synchronize their beat by means of hydro-mechanical 
coupling. 

The hydrodynamic flows generated by beating flagella are characterized by low 
Reynolds numbers [8 61 6^. The Reynolds number characterizes the relative impor¬ 
tance of inertial over viscous forces, and is defined as Re = pvL/rj given a typical 
speed V and size L of active motion. Here, p and p denote the density and the dy¬ 
namic viscosity of the fluid. In the limit of zero Reynolds number, one obtains the 
Stokes equation, which for an incompressible, Newtonian fluid reads 


63 


0 = —Vp + pV'^v. 


(18) 


Here, p and v denote pressure and flow velocity. For the following, it will be im¬ 
portant that the Stokes equation is linear and time-reversible, i.e. the fluid flow v 
will simply change its sign if a flagellum would play its swimming stroke backwards 
in time. This time-reversibility has well-known consequences for self-propulsion of 
microswimmers at low Reynolds numbers, implying that a microswimmer must use 
a non-reciprocal swimmi ng stro ke that explicitly breaks time-reversal symmetry to 
allow for net propulsion K 62 . Similar arguments also apply to synchronization by 


hydrodynamic forces in the Stokes limit. 


3.1 Lack of hydrodynamic synchronization in a minimal model 


Minimal models that retain key symmetries of flagellar swimming and synchroniza¬ 
tion have proven very useful to elucidate fundamental principles of hydrodynamic 
synchronization 38 47 48 51 64 -66 . Specifically, a beating flagellum has been mim¬ 
icked by a sphere moving along a perfect circle 47 , see Figure]^. This description is 
motivated by the observation that each point on a flagellum traces out a circular orbit 
with respect to the reference frame of the cell while the flagellum beats. An entire 
flagellum could be described as a collection of several of such idealized spheres 


30 


We consider the case of two spheres of radius a that are constrained to move 
along respective circular trajectories of radius A, ri((pz,) and which are 

parametrized by phase angles pL and pji and separated by a distance d, see Figure 


rn(pH) 


A 


cos PR 
sin pr 


(19) 


Each sphere is driven by a driving force Fj of constant magnitude /o that acts tan¬ 
gential to its circular track, i.e. Fj = /ot^-, for j e {L,i?}, where tj = r~^dTjldpj 
denote the local tangent vectors. In the inertia-less limit, these driving forces are 
balanced by the hydrodynamic friction forces Hl and H/j exerted by the spheres on 
the surrounding fluid. 


Fi-tL=Hi-tL, Fij • = Hi^ • tij. (20) 

Note that force balance between driving force and hydrodynamic friction force holds 
only tangential to the track. Along the normal direction, there can be a difference of 
forces. In fact this force difference amounts to the constraining forces needed to keep 
the spheres on their circular track. We use short-hand notation Pj = AHj ■ tj and 
Qj — AFj ■ tj for j e {L, R}. 

In the limit of small spheres, a <C A, and large separation distances, d ^ A, the 
hydrodynamic friction forces are readily computed using the Stokes drag coefficient 
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Fig. 3. Lack of hydrodynamic synchronization in a minimal model. A. Each point on a 
beating flagellum traces out a circular orbit. To study qualitative aspects of hydrodynamic 
synchronization, we can replace each beating flagellum by a sphere that moves along a 
circle. B. Two co-rotating spheres of radius a move along circular tracks of radius A, rr and 
r_R, which are separated by a distance d. Active driving forces of constant magnitude act 
on the spheres tangential to their track. Additionally, constraining forces keep the spheres 
on track. C. The motion of the two spheres induces flow of the surrounding fluid (here 
indicated in color-code). Thereby, the motion of the two spheres becomes coupled by virtue 
of hydrodynamic interactions: the flow induced by the motion of one sphere propagates to 
the second sphere, where it exerts a hydrodynamic friction force. D. Although hydrodynamic 
interactions induce a weak coupling between the instantaneous motion of the two spheres, 
the net effect of this coupling is zero. This can be seen from symmetries: a reflection Mx 
at the x-axis, and a time-reversal T both change the dynamics of the spheres in the same 
way. As a time-reversal flips the stability behavior of synchronized states, while a reflection 
leaves this unaltered, we conclude that synchronization cannot occur in this minimal model, 
see text further explanation. 


A 




o g 




B 


^ c* 

fM, 

o o g 



Fig. 4. The minimal model revisited. The spatio-temporal symmetry of the minimal model 
of two co-rotating spheres has to be broken to allow for net synchronization by direct hy¬ 
drodynamic interactions. A. A stationary third sphere (gray) breaks the dynamic equiva¬ 
lence between the system reflected by Mix and the system after time-reversal T. However, 
a reflection Mly at the j/-axis again yields a system that is dynamically equivalent to the 
time-reversed system. This second spatio-temporal symmetry rules out net synchronization 
for a pair of co-rotating spheres, if the third stationary sphere is in a symmetric position 
as drawn here. B. The situation is different for the case of counter-rotating spheres, which 
more closely mimics the mirror-symmetric beat of the two flagella of the Chlamydomonas 
cell shown in Figure]^ In this case, either mirror operation yields a dynamics that is differ¬ 
ent from the time-reversed dynamics. Indeed, synchronization by hydrodynamic interactions 
does occur [^. 
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7 = dnija of a single sphere, and the Oseen tensor G(r) that describes how the flow 
generated by the motion of one sphere propagates to the other sphere, 


'-+rirj\r 


( 21 ) 


Hl = 7 [tcl - G(ri^ - r^) • {AipF.tR/-i)] , 

and similarly for H^. The Oseen tensor is given by Gij(r) = (87r?7)“^(| 

For the following, all we need is that the generalized friction forces Pj can be written 
as a linear function of the phase speeds pj, e.g. Pl = Pll/^l + PlrP’r- This general 
form holds true without the approximation of small spheres, for more general geome¬ 
tries, and even generalizes to free-moving swimmers 
balance equation, Equation (20), as 


50 . Thus, we can write the force 


T'{ipL,‘PR) 


pR 


Ql 

Qr 


( 22 ) 


This force balance equation describes two coupled oscillators and can be re-cast in a 
form similar to Equation § studied above. 

We can understand the synchronization dynamics solely by considering symme- 
let T be the time-reversal operator 


tries 49 


while M-a; and M.y denote spatial 


mirror operations at the x- and y-axis, respectively, see Figure Since there is no 
explicit time-dependence in the Stokes equation. Equation (ll^, and the laws of hy¬ 


drodynamics are invariant under reflections, the friction matrix T{ipF,ipR) must be 
invariant under these symmetry operations. These symmetry operations change the 
phase speeds pj and generalized driving forces Qj. 

The symmetry operations T and induce the same transformation of the dy¬ 
namic equations of motion, Equation \22\ . Initial conditions are transformed differ¬ 
ently, though, Mx ■ i^j(O) —>■ —ipj{t)) and T : <^5^(0) —>■ ‘/^^(O), yet this does not change 
the argument. Thus, applying either Aix or T, we should find the same stability be¬ 
havior of synchronized steady-states. The crucial point is now that a time-reversal T 
will flip the stability behavior of any synchronized steady-state: a stable steady-state 
will become stable and vice versa. In the Adler equation. Equation ([^, this corre¬ 
sponds to A —> —A. On the other hand, a spatial mirror operation will not change the 
stability behavior of any steady-state. We conclude that there cannot be any stable 
or unstable steady-state of the untransformed system. In fact, all steady-states are 
neutrally stable, corresponding to the case A = 0 in Equation 

Thus, the symmetries of the Stokes equation rule out stal^ synchronization in 
this minimal model. A number of specific generalizations have been proposed, all of 
which break the TAf ^-symmetry in a specific manner as discussed in the next section. 


Note on generalized forces. The quantities Pj above represent generalized forces 
with units of a torque, which are conjugate to ipj in the sense of Lagrangian mechan¬ 
ics of dissipative systems [M 67 ^ . Formally, these generalized friction forces are 
defined as Pj — dTZ/dpj, where TZ denotes the total rate of hydrodynamic dissipa¬ 
tion, which plays the role of a Rayleigh dissipation function. Lagrangian mechanics is 
particularly useful to describe mechanical systems with constraints and applies both 
in conservative and non-conservative systems with quadratic dissipation function. 


3.2 Hydrodynamic synchronization by broken symmetries 


Already Taylor suspected that flagellar synchronization is a hydrodynamic effect 
yet the first experimental demonstrations were only recent 


30 


39] . Again, simple 

model systems served as a proof of principle that synchronization by hydrodynamic 


38 
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forces is possible, both in theory ^ 48 51 64 -66 and experiment 42 -46 . These 


model systems demonstrated equally clearly the need for non-reversible phase dy¬ 
namics and thus broken symmetries 49 . In these models, symmetries were broken 
either by wall effects 47 , by additional elastic degrees of freedom 48 64 , or by 


evoking phase-dependent driving forces 44,51 


Symmetry breaking by boundary walls. In the original formulation of the two- 
sphere-model by Vilfan et al. 47 , a no-slip boundary close to the two circular trajec¬ 
tories was introduced, which changes the friction matrix F, and in particular breaks 
the symmetry of F with respect to the spatial mirror-operation Aix- Later work used 
a third, stationary sphere, instead of the boundary wall Note that in addi¬ 


tion to the TAl a;-symmetry, one can make a similar symmetry argument also for a 
reflection M.y at the y-axis. To break also this second symmetry, elliptic trajectories 
with a tilted orientation with respect to the boundary wall were considered in the 
work cited to provide stable synchronization 47 . Interestingly, the second symmetry 


argument only applies to the case of co-rotating spheres, but not for counter-rotating 
spheres with lol = —ojr, which makes synchronization easier in that case 40 50 , see 
Fig.|4] 


Symmetry breaking by amplitude compliance. Lenz et al. introduced an elastic 
compliance into the model by considering the radii AL{t) and of the left and 

right track as additional degrees of freedom. Specifically, the spheres were held by 
elastic springs with rest length A and spring constant k, which can slightly change 
their length in response to hydrodynamic forces. The synchronization strength was 
found to be inversely proportional to the stiffness of the springs, A ~ 1/k. This 
result is consistent with zero synchronization strength A —>■ 0 in the limit of a hard 
constraint A = Aq for fc —> oo. Several experimental realizations highlighted the 
role of elastic compliance for hydrodynamic synchronization in pairs of man-made 
oscillators, includin g he lices rotating at low Reynolds numbers in a tank of highly 
viscous silicone oil 69 , or micro-rotors driven by laser light 45 . A recent study 


of the synchronization of beating flagella by direct hydrodynamic interactions also 
suggested an important role of elastic waveform compliance of the flagellar beat for 


synchronization 39 


Symmetry breaking by phase-dependent driving forces. Golestanian et al. dis¬ 
cussed phase-dependent driving forces Qj{^j) as a generic mechanism for hydrody¬ 
namic synchronization. Such phase-dependence represent another way to break TM.x- 
symmetry. In computational models of flagellar synchronization for the bi-flagellated 
green alga Chlamydomonas, phase-dependent driving forces were indeed found to 

Sim 


contribute to synchronization 


In short, synchronization arises from a “reso¬ 


nance” between a phase-dependence of the active driving and the phase-dependence 
of the hydrodynamic friction coefficients. In a coarse-grained theoretical description 
of flagellar beating, phase-dependent driving forces of the flagellar beat can be com¬ 
puted from the phase-dependent rate of hydrodynamic dissipation 30 m. Figure 


shows the phase-dependence of the hydrodynamic dissipation rate computed for a 
free-swimming Chlamydomonas cell. 


Symmetry breaking by inertial effects. 

but non-zero Reynolds number Re 


52 


Theers et al. considered the effect of small, 
The presence of unsteady acceleration terms 
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Fig. 5. Phase-dependent hydrodynamic dissipation. The rate of hydrodynamic dissipation 
for a swimming Chlamydomonas cell changes during the flagellar beat cycle, suggesting 
phase-dependent flagellar driving forces. The green curve shows the total hydrodynamic dis¬ 
sipation, while the blue curve depicts a suitably defined dissipation in the far field [70| . The 
insets in the lower row show the flagellar shapes at the different phase of the beat cycle 
together with a numeric computation of the flow profile. The total hydrodynamic dissipa¬ 
tion rate is maximal during the recovery stroke, when the two flagella move close to the cell 
body. In contrast a hydrodynamic dissipation rate computed for the far-held of huid how 
is maximal during the effective stroke. In a simple description of the hagellar beat, the hy¬ 
drodynamic dissipation rate corresponds to an active driving force of the hagellar beat [30| . 
Phase-dependent active driving forces represent a general mechanism for hydrodynamic syn¬ 
chronization that explicitly breaks symmetries that rule out synchronization 
from 


70 


51 . Modihed 


in the Navier-Stokes equation breaks it time-reversal symmetry and thus allows for 
stable synchronization. The synchronization strength scales as A ^ (Re)^/^. A similar 
line of argument should also apply to complex fluids with visco-elastic properties or 
nonlinear constitutive relations, which likewise break time-reversal symmetry. 


3.3 Flagellar synchronization independent of hydrodynamic interactions 


In addition to mechanisms for flagellar synchronization that rely on direct hydrody¬ 
namic interactions between two flagella as discussed in the previous section, Friedrich 
et al. recently introduced an alternative mechanism that is independent of hydro- 
dynamic interactio ns, but relies instead on a coupling between flagellar beating and 
swimming motion 30 40 . It was proposed that this synchronization mechanism ap¬ 


plies in particular to free-swimming Chlamydomonas cells. In short, both the left 
and right flagellum of a Chlamydomonas cell were modeled as phase-oscillators with 
respective phase, (pL and (pR. The speed of the beat ipj, j G {L,i?}, depends on the 
instantaneous hydrodynamic load acting on the flagella, which in turn depends on 
the swimming motion of the cell, characterized by its rate of rotation d 


(pL =uJo- Kv’L)a, 

(pR = Wo fJ,{(pR)d. 


(23) 

(24) 


These dynamical equations can be derived from simple assumptions using a general 
framework of force balances. In particular, the coupling function fj, can be computed 
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for any given beat pattern in terms of hydrodynamic forces acting on the two flag¬ 
ella, see Figure Direct hydrodynamic interactions between the two flagella can be 
accounted for, but have been neglected in Equation (23) for sake of simplicity. Asyn¬ 


chronous beating of the two flagella causes a rotation of the cell with rotation rate 
a 

pi^PL, PR)a = v{ipl)pl - v{ip)rlpr. (25) 


Equation (25) represents a torque balance equation that balances the torques for a 


rigid body rotation of the entire cell (on the left-hand side) and the sum of the torques 
generated by the beat of the left and right flagellum, respectively (right-hand side). 
Note that Equation 25 implies d = 0 if both flagella beat in-phase, as expected by 
symmetry. 


We can re-arrange the dynamical system given by Equations (23) and (25) into 


the form of an Adler equation for the phase difference 5 = (pR — Pr, see Equation 
This allows to express the effective synchronization strength A in terms of the various 
coupling functions 30 


A = - 


dp 


2fi{py'{p) 


pip,p) - 2fi{(p)iy((p) 


(26) 


The coupling functions v, and p can in turn be computed in terms of hydrodynamic 
friction forces. For wild-type beat patterns, it was calculated A = 0.3, which implies 
stable in-phase synchronization 30 . Note that alternative beat patterns can result 
in A < 0, which renders the in-phase synchronized state unstable, and favors instead 
anti-phase synchronization with i5 « tt. 


3.4 The beating flagellum is a noisy oscillator 


The flagellar beat represents a biological oscillator with small, yet perceivable fluc- 

By mapping the periodic shape deformations of the 


tuations 


27-29,34,57,71,72 


flagellar beat on the normal form of a Hopf oscillator. Equation (1^, it is possible to 


quantify both phase and amplitude fluctuations 57 . In the simplihed case of a noisy 
phase oscillator, given by Equation ([^, we can compute the quality factor of noisy 
oscillations 

The quality factor is a measure after how many oscillation cycles the oscillator phase 
decorrelates due to noise. For bull sperm flagella, it was found Q = 38.0± 16.7 (mean± 
s.e.) [57| . This value is on the same order-of-magnitude as previous, indirect measure¬ 
ments for Chlamydomonas flagella based on the frequency of phase-slips in pairs of 

The corresponding noise strength D^p = u}q/{2Q) is several 


synchronized flagella 28 


orders of magnitude larger than an estimate for contribution from passive, thermal 
fluctuations based on a fluctuation-dissipation theorem 57 . Thus, fluctuations of the 


flagellar beat are of active origin Fast fluctuations of the flagellar beat with 

correlation times on a millisecond time-scale have been attributed to small-number- 
fluctuations in the collective dynamics of the ap proximately 10^ molecular motors 
inside the sperm flagellum that drive its beat [^. Thus, fluctuations of the flagellar 
beat represent mesoscopic signatures of stochastic activity of individual molecular 
motors. Subsequent, more refined analyses of flagellar fluctuations have shown that 
amplitude fluctuations are phase-dependent, being minimal during the initiation of 
bends at the proximal end of the flagellum 34,72 . Additionally, slow fluctuations of 


the flagellar beat with correlation times on a second time-scale have been observed, 
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Fig. 6. Dynamic force-velocity-relation of the flagellar beat A. In a theoretical description 
of flagellar swimming and synchronization in the green alga Chlamydomonas, each beating 
flagellum was described by a phase oscillator with respective phases ipr and ipn [30| . In this 
description, the flagellar phase speed changes in response to a change in hydrodynamic 
load, which can in particular result from a rotation of the entire cell with rotation rate a. 
In the absence of any rotation, d = 0, the flagellum beats at its intrinsic frequency ljJl- For 
a rotating cell, additional hydrodynamic friction forces act on the flagellum to which the 
flagellum responds by a change in phase speed with susceptibility fj,{g>). This susceptibility 
characterizes a dynamic force-velocity relationship of the beating flagellum. B. The coupling 
coefficient can be computed in terms of hydrodynamic friction forces, assuming that the 
flagellar oscillator operates at full load with negligible internal dissipation. The susceptibility 
has been measured in experiments with Chlamydomonas cells that displayed oscillatory 
rotational motion due to a lack of flagellar synchronization. During the effective stroke of 
the left flagellum (93 ~ 0 ), a counter-clockwise rotation of the entire cell (d > 0 ) increases 
the motion of the flagellum relative to the fluid at rest. As a consequence, the hydrodynamic 
friction forces acting on this flagellum increase, and the flagellum slows down (/r > 0 ), both 
in theory and experiment. During the recovery stroke, the effect is opposite: The active 
shape change of the flagellum and a passive counter-clockwise rotation partially cancel each 
other. This reduces the hydrodynamic friction forces acting on the left flagellum, causing it 
to speed its beat up {fi < 0). Modihed from 


30 


which may reflect noise in chemical signaling cascade that control the shape of the 
flagellar beat 


34 


The synchronization of several flagella represents a stochastic phenomenon that 
arises from a competition between a mechanical coupling (which favors phase-locking), 
and the effects of noise (which counter-act synchronization) [27f 29 57 . 


4 Discussion 


Flagellar synchronization represents a model system for the collective dynamics of 
biological oscillators. Experiments have demonstrated that a mechanical coupling 
by viscous forces can indeed synchronize several flagella 

each of 


Theory has outlined 

different physical mechanisms for hydrodynamic synchronization 47|48 51|52 
which breaks the time-reversal symmetry of hydrodynamics at low Reynolds r 


symmetry of hydrodynamics at low Reynolds numbers 
in a different way. Which of the different mechanisms dominates in different biological 
systems will have to determined by future research. Open questions relate in particular 
to the role of elastic anchorage of the flagellar apparatus 37 , and the waveform 
compliance of the flagellar beat 34 . It is not known, which features of the flagellar 


beat patterns determine whether in-phase or anti-phase synchronization will be stable, 
as observed e.g. in flagellar mutants 


31 


Previous theory suggested that states of synchronized dynamics correspond to 
either a maximum or a minimum of the rate of hydrodynamic dissipation [49| . How¬ 
ever, it is not known if this rule is universally applicable and which features determine 
which extremum corresponds to a stable synchronized state. 
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We are only beginning to understand how the flagellar beat responds to time- 


varying external forces 30 34 . Such an understanding will be important not only 


to understand synchronization in collections of beating flagella by mechanical cou¬ 
pling, but may also be informative o n the m icroscopic mechanisms of motor control 
that regulate flagellar bending waves 13 -15 . Future research can bridge the gap be¬ 


tween educative minimal models of flage. 
biological systems. 


far synchronization, and the complexity of 
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